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Abstract 

Based on recent results obtained by the authors on the inverse scattering 
method of the vector nonhnear Schrodinger equation with integrable boundary con- 
ditions, we discuss the factorization of the interactions of A^-sohton solutions on the 
half-line. Using dressing transformations combined with a mirror image technique, 
factorization of soliton-soliton and soliton-boundary interactions is proved. We 
discover a new object, which we call reflection map, that satisfies a set-theoretical 
reflection equation which we also introduce. Two classes of solutions for the re- 
flection map are constructed. Finally, basic aspects of the theory of set-theoretical 
reflection equations are introduced. 
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1 Introduction 

In this article, we consider the Manakov model fT] or more generally the vector nonlinear 
Schrodinger equation (VNLS) for the n-component complex- valued vector field R{x,t), 
in the focusing case, 

,_^_^2B.B^R^,, (1.1) 

(1.2) 

where t) is the transpose conjugate of R{x,t), and we restrict it to the half-line 

X > hj imposing integrable boundary conditions derived in [2j: Robin boundary 
conditions 



R^{0,t) -2aR{0,t) = , aeR, (1.3) 

or a mixture or Neumann and Dirichlet boundary conditions 

R,{0,t)=0 , jeS , (1.4) 
i?fcx(0,t)=0 , ke{l,...,n}\S, (1.5) 

where S" is a subset of {1, . . . ,n}. Supplemented with the initial condition R{x,0) = 
Ro{x), we therefore consider an integrable initial-boundary value problem^. In [2j, the 
inverse scattering method valid for the problem on the full line has been suitably adapted 
to yield reconstruction formulae for the function R{x, t). In particular, explicit A^-soliton 
solutions have been obtained. An interesting effect has been identified when a soliton 
bounces off the boundary: its polarization can be altered and a redistribution of the 
amplitudes of the various components occurs through a process of transmission between 
modes. However, a full understanding of the dynamics of the solitons on the half-line is 
still lacking and it is the object of this paper to discuss this point. 

For the problem on the line, this question has only been studied rather recently in 
[H|5]. The fundamental result is that soliton interactions (or collisions) factorize in the 
vector case like in the scalar case. In the scalar case, factorization of soliton interactions 
is a well-known and important fact. Roughly speaking, it means that as time evolves 
from — oo to oo, two solitons with given (different) velocities and amplitudes as t — )• — oo 
will collide and recover their initial properties (velocities, amplitudes, shapes), the effect 
of the interaction being only position and phase shifts. This fact generalizes to iV-soliton 
collisions. In this case, the resulting velocities and amplitudes of the solitons as t — > cxo 
is independent of the way the collisions took places. This gives rise to the notion of 
factorization of soliton interactions. In the vector case, solitons are specified by their 
so-called polarization vectors on top of their velocities and amplitudes, and the property 
of factorization now involves these polarizations. It is therefore a highly nontrivial fact 
that factorization still holds. In terms of the polarizations, this means that the set 

^Of course, initial and boundary values are supposed to be compatible and, to fix ideas, we work in 
the space of Schwarz functions which is known to be appropriate for the inverse scattering method, see 
e.g. t3] 



1 



of polarization vectors of the solitons obtained as t — ^ oo is independent of the way 
the colhsions between sohtons took place but depends only on the initial polarization 
vectors. 

This property in the vector case is intimately related to the existence of a so-called 
Yang-Baxter map which satisfies the set-theoretical Yang-Baxter equation (YBE) as 
originally proposed in [6j and studied in more details in jTllBllS]- This map describes the 
relation between two polarization vectors before and after the collision. The proof of the 
factorization property was established using two related but rather different routes in |1] 
and [5]. In j3], the final result of the inverse scattering method, i.e. the explicit A^-soliton 
solution, is carefully studied in the limits t — )■ ±00 to extract the factorization property 
by looking directly at the evolution of the polarization vectors. In [5], the results obtained 
originally by Manakov by asymptotic analysis of the inverse scattering method are used 
to construct an inductive argument on the number of solitons. Then, this is cast in the 
formalism of Yang-Baxter maps and the associated Lax pairs to confirm factorization 
using matrix polynomials refactorization properties. The Yang-Baxter map is realised 
in this context by a map acting on pairs of polarization vectors. It is the mathematical 
representation of a soliton-soliton interaction. 

In essence, the set-theoretical YBE ensures the factorization property. The situation 
is analogous to the quantum case where factorization is ensured by the quantum YBE 
but we emphasize that the nature of the objects that are involved is fundamentally 
different. In fact, the set-theoretical YBE involves maps on cartesian products of sets 
and provides a very general framework. The quantum YBE can be seen as a special case 
of general theory. 

By studying the question of factorization in the presence of a boundary, we discover 
a new framework which, jointly with the set-theoretical YBE, ensures the factorization 
property on the half-line: the set-theoretical reflection equation. We also find two classes 
of solutions for this equation which we call reflection maps. 

The paper is organised as follows. In section El we use the formalism of the Riemann- 
Hilbert problem and the dressing method [TOl [TT] to establish Theorem 12.71 which is 
the central technical tool in our approach. First, it allows us to rederive the above 
factorization property on the line directly at the level of the dressing method. On top 
of providing a fresh look and a new proof of the result, this allows us to introduce all 
the relevant definitions, notations and results needed for our purposes. Then, in section 
[3l we combine these results with the mirror image technique developed in [2] to prove 
the factorization of the soliton-soliton and soliton-boundary interactions. This gives 
rise to the introduction of a new object, the reflection map, which is a map acting on 
one polarization vector at a time and is the mathematical representation of the soliton- 
boundary interaction. Just like the factorization of soliton-soliton interactions is related 
to the set-theoretical YBE, we find that the factorization property with a boundary is 
related to an equation which we call the set-theoretical reflection equation. This and the 
construction of explicit reflection maps form the main result of this paper presented in 
Theorem 13.31 Section H] is then devoted to the presentation of a few basic elements of 
the theory of set-theoretical RE in an abstract setting. A few technical arguments are 
collected in appendices. 

It is remarkable that our results complete the formal analogy that exists between 
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quantum integrable systems and classical soliton partial differential equations. We al- 
ready mentioned this analogy at the level of the YBE. Here, our findings provide a 
further link with the reflection equation that appears in the theory of quantum inte- 
grable systems with boundaries [T2| [13]. In the particular case of VNLS, the complete 
classification of the solutions of the quantum reflection equation was obtained in [14j . 



2 Dressing method, factorization and Yang-Baxter maps 

In this section, we first review some basic concepts of the inverse scattering method 
(ISM) and its relation to the dressing method. The main contribution relies on the 
permutability property of dressing transformations which is stated in Theorem 12.71 It 
then allows us to provide a new proof of the factorization property for VNLS on the line 
and to show the underlying Yang-Baxter structure. 

We refer readers to e.g. j3l [15l [16] for more detailed presentations of the ISM, and 
to e.g. [T7] for a presentation of the dressing method. 



2.1 Lax pair formulation 

Define Q{x,t) as the following [n + 1) x [n + 1) matrix-valued field 

The VNLS equation (11. II) is the compatibility condition ($^.f = $tx) of the following 
linear problems for the matrix-valued function ^{x,t, 

$, + ^A;[S3,$] = Q$, (2.2) 

$j + 2«A;'[S3,<l>] = gr$, (2.3) 

where 

^3 = ( o" , Qt = 2kQ - iQ, S3 - tQ^ S3 , (2.4) 

/„ being the n x n identity matrix. Equations (12. 2[ 12. 3p form the Lax pair formulation 
for VNLS. One observes that Q satisfies 

g=-Q^ (2.5) 

where ''' denote the transpose conjugate operation. This implies that, ^{x, t, k) being so- 
lution of the Lax pair, $''"(x, t, k*) satisfies the same equations as "^{x, t, k) = $~^(x, t, k) 
i.e. 




^, + zA;[S3,^] = -^Q, (2.6) 
^t + 22p[S3,^] = -^Qr- (2.7) 



^From now on, we drop the x, t and k dependence for conciseness unless there is ambiguity. 
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We define two Jost solutions X and Y of tlie Lax pair (12. 2^ 12. 3p satisfying 

lim e''^(^'*''=) ^'X{x, t, A;)e-^'^(^'*'^) =J„+i , keR, (2.8) 

X— — oo 

lim e^'^^^'*''^) t, fc)e-^'^(^'*'^) =J„+i , A; e M , (2.9) 

x—^oo 

where (j){x, t, k) = kx + 2k'^t. They enjoy the following properties 

det X{x, t, k) = det y (x, t,k) = l. (2.10) 

• X and Y can be splitted into the following "column" vector^ form 

X = {X+,X-), Y = {Y-,Y+), (2.11) 

where y+ (resp. X~, Y^) are analytic and bounded in the upper (resp. 
lower) half fc-complex plane, which are denoted by ( resp. C^). 

• X and Y are related by the so-called scattering matrix S{k) 

X{x, t, k) = Y{x, t, fc)e-*<^(^'*''=)^3 S{k) e*<^(^'*''=)^^ , keR, (2.12) 
where S{k) can be splitted into block matrices of natural size^ 



b+{k) a-{k) 

where a^(/c) allows for analytic continuation into respectively. 

X-\x,t,k) = X^{x,t,k*) , Y-\x,t,k) = Y\x,t,k*), (2.14) 

and hence 

S{k)-^ = S^ik*). (2.15) 

In components, we denote 

Accordingly, c^{k) has an analytic continuation into C^. In particular, one has 
the following relations 

(a±)t(r) =c^(A;), keC^, (2.17) 
deta+(A;) = c+(A;), k e C+ . (2.18) 



•^Here, the left "column" vector is made of the first left n columns and the right one is made of the 
remaining column. This column vector representation will always be used in the rest of this paper. 
^For instance, a'^ is an n x n matrix while a~ is a scalar. 
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For VNLS, there are two equivalent sets of scattering functions {a^ (k) , (k)} and 
{c^ (k) , (k)} available to reconstruct R{x,t) in the inverse part of the ISM [16|. In 
the rest of this paper, we choose to work with a^(/c) and b^{k). The information on the 
soliton solutions depends on the analytic structure of a^{k). 

A key observation in the development of soliton theory is that the ISM can be cast 
into a Riemann-Hilbert problem. The latter offers a natural and powerful framework 
for the so-called dressing method [T0| ITT] . In the next subsection, we present general 
results about the dressing method before returning to the case of VNLS. 

2.2 Riemann-Hilbert problem and dressing method 

We mainly follow the presentation of chapter 3 of [T7] and begin by stating the main 
propositions about the connection between Riemann-Hilbert (RH) problems with zeroes 
and the dressing method. We refer readers to fTT] for details and in particular for the 
proofs of Propositions 12.31 and 12.41 below. 

Consider the following matrix RH problem with canonical normalization 

J^{k)J-{k) = J{k) , keR, Urn J^{k)^I, (2.19) 

)-oo 

where J^{k) is the jump matrix satisfying detj'{k) ^ for /c G M. This RH problem 
is discussed in some details for instance in section 7.5 of [TH]. The matrix J^{k) is 
analytic in C^. This problem has a unique regular solution J^ik) that is a solution 
with det J^{k) 7^ 0. In the construction of soliton solutions, the notion of RH problems 
with zeroes plays an important role. Although not the most general one (see e.g. |19]). 
the following definition is sufficient for our purposes. 

Definition 2.1 A matrix M{k) is said singular at k = ko if det M{kQ) = and if in 
the neighbourhood of k^ 

M{k) = Mo + {k~ ko)Mi + 0{k - kof , M-\k) = + N^ + Oik - ko) .(2.20) 

Definition 2.2 A RH problem with zeroes at k^ G C"*^, j = 1, . . . , N is a RH problem 
as in fl2.19p where J'^{k) is singular at k^, j = 1, . . . , N . 

Then one proves 

Proposition 2.3 Fixing the subspaces Vj = Im J'~^{k)\,,+ andUj = Ker {k)\,^,- , 
j = 1, . . . , N determines uniquely the solution of the RH problem with zeroes at k^ G C^. 

In general, there is no known closed-form formula to solve a matrix RH problem. 
However, once a regular solution is known, it is possible to construct singular solutions 
from it. 

Proposition 2.4 Let Jl^{k) be the singular solution at k^ G C"*^ with 

lmJ+{k)\^_^^ = Vo, KeTj-ik)\^_^,=Uo, (2.21) 
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and let J^{k) he the solution of the same RH problem regular at k^. Then J^{k) can 
he written as 

J^{k) = J-+(A;)(J„+i + (^-I^)no) , J~{k) = (I^^, + (^fl^)no)J,{k) . (2.22) 

Here Uq is the projector defined by 

KerHo = (Jo^(A:o+)) Vo , ImHo = Jo'{ko)Ko • (2.23) 

Proposition 12.41 introduces what is called a dressing factor (of degree 1) which transforms 
J^{k) regular at into J^{k) singular at k^. This gives an algorithm to construct 
a singular solution J^{k) at distinct k^ G C^, j = 1, . . . ,N from a regular solution 
J^{k). Precisely, let k^ G C^, j = 1, . . . , and the corresponding subspaces Vj, Uj 
be given. Use Proposition 12.41 repeatedly to construct J^{k) singular at k"^ recursively 
from J^{k) starting from A;f , A;^ up to k^. Consequently, J^{k) can be written as 

j^{k) =X{k) + ^^ni) . . . + ^^n^) , (2.24) 



J-{k) = I + ^^n^) . . . f/n+i + ^^ni ) j,-{k) , (2.25) 

TV / V 1 



where, for j = 1, . . . , 



(2.26) 



imn, = ( + n.-i j • • • [in+i + ^^^i^ni j Jo{kj)u, . (2.27) 

Now comes a simple but fundamental observation which is absent in chapter 3 of |17) . 
In the above construction, one can iterate Proposition 12.41 by using a different order on 
the k^. Let Sn be the permutation group on the set {1, . . . , N} and let cr G Sn- Denote 
the image of (1, . . . , A^) under a by (o"(l), . . . , 0"(iV)) and introduce = k^^^-y Then, 
the subspaces corresponding to are Vq-q), tia(j)- Repeating the previous procedure, 
starting from up to k^, we obtain 

'-^^'n-V.. f/^,, + 7 '-^"^ n^ 1 , (2.28) 
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where, for j = 1, . . . , iV 



(2.30) 

-0) -0-1) / V -u) -(1) / 

(2.31) 



It can be checked by direct calculation that Vj = Im J^{k) and Uj = Ker J' (k) 



k=k] 



k=k- 



j = 1, . . . ,N SO that Proposition 12.31 implies that J^{k) = J^{k). In turn, this implies 
that the product of dressing factors in (I2.28p . f l2.29p is equal to the product of dressing 
factors in f l2.24p . (I2.25p . This construction introduces the notion of dressing factor of 
degree which transforms J^{k) into J^{k) singular at /cj*^, j = 1, . . . , A^. Proposition 
12.31 also ensures that a dressing factor of order factorises into dressing factors of 
degree 1 and that the order of the factorization is irrelevant. At this stage, we note that 
this fact is known in various disguises (e.g. as the Bianchi permutativity property) and 
in various context (e.g. in the theory of matrix polynomials) but, to the best of our 
knowledge, it has not been presented anywhere in the above fashion. This is why we 
formulate it as a theorem below in a form suitable for our purposes. 

Remark 2.5 It is important to realize that this does not mean at all that the individual 
factors in a dressing factor of degree N commute. Indeed, in general IIJ ^ '^cr{j) ■ This 
can happen in special circumstances and in that case, the interaction of the solitons is 
trivial from the polarization point of view. The important message here is that, in the 
factorization of a dressing factor of degree N , the equations governing the projectors are 
crucial, in particular the order in which they appear is important. With this in mind, 
we introduce a notation that will help us formulate the main theorem of this section and 
prove the factorization of soliton interactions in the following sections. 

Definition 2.6 Given J^{k) a regular solution of the RH problem fl2.19p . Let a G Sn 
he given and write (cr(l), . . . ,a{N)) = (ii, . . . ^in)- Given and Vj, Uj, j = 1, . . . , N , 
a general dressing factor of degree 1 is defined as, for I < i < N , 

A,,{ii...i,_i}(/i;) = in+i + j^J n^,,fa...»,_i} (2.32) 

where 

Kern,„{,,..,,_,} = [A,(/c+)...A,_„{.,.....})(A:+)]"' W(A:+)) V., , (2.33) 

Imn,„{,,..,,_,} = [A.(A:,:;)... Jo"(A;,^)W.,. (2.34) 

Finally, the dressing factor of degree N is denoted as Di j^[k). 
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Note that in the case i = 1, we denote |}(A;) = Di^{k). This convention will be 
adopted in the rest of the paper for any quantity involving sets of indices as subscripts. 
With this definition and the understanding from the above discussion, we have proved 
the following 

Theorem 2.7 A dressing factor of degree N can he decomposed into N\ equivalent prod- 
ucts of N dressing factors of degree 1 

D^,„N{k) = D,,{k) . . . A^,{n....^^i}(A;) , (2.35) 

where (zi, . . . , i^) is an arbitrary permutation of {1, ... , N). 

In the next subsection we want to use the results of this subsection for a RH problem 
arising from a Lax pair formulation. The latter involves the variables x and t so that in 
fact, we are dealing with a parameter dependent RH problem of the form 

J^{x,t,k)J-{x,t,k) = J{x,t,k) , keR, Mm J{x,t, k) ^ I . (2.36) 

The success of the dressing method is related to the fact that all the results seen in this 
subsection go through for this parameter-dependent RH problem provided one works 
with (x, t)-dependent subspaces Vj{x, t), Uj{x, t), j = 1, . . . , N (and the associated pro- 
jectors) which are simply related to Vj, Uj, j = 1, . . . , N, by 

Vj{x, t) = ^{x, t, k+)Vj , Uj{x, t) = ^{x, t, k-)Uj , (2.37) 

where if is a solution of the so-called undressed Lax pair equation^. 



2.3 Application to VNLS: reduction 

We begin by collecting some known facts that can be gathered for instance from jSl [IHl 
|2U] . The scattering system defined in fl2.12[) can be rewritten as the following {x,t)- 
dependent RH problem 

J^{x,t,k)J'{x,t,k) = J{x,t,k) , keR, Mm J^{x,t, k) ^ In+i . (2.38) 

|A;|^oo 

This is achieved by defining 

J+(x, t, k) = (^^^^^^ ((X+, r+) ' (x, t, k) , J-(x, t, k) = {Y-, X-) (x, t, k) , 

(2.39) 

and 

J(x, t, k) = e-*'^("'*''=)^3 (^^/^^^ ei0(-Afc)S3 ^ ^ g K _ (2.40) 

In particular, 

det J+(x,t,fc) = deta+(fc) , det J"(x, t. A;) = a-(A;) . (2.41) 

One then has 

^Indeed, (p satisfies Uip = (fx, Vip = (ft, with Ut — Vx + [U,V] — (zero curvature condition). The 
Lax pair (j2.2l 12. 3p is another way to express the zero curvature condition. For VNLS, if is simply 
related to $ which satisfies (E^IOl) hy (p = $e-»</'(^^*>fe)S3, 
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Proposition 2.8 Consider the RH problem defined in (I2.38p . Let J^{x,t,k) be a so- 
lution of the Riemann-Hilbert problem, then J~^{x,t,k) (resp. J~{x,t,k) ) satisfies the 
Lax pair \2.(X \2. 7|) (resp. 112. 2. \2.3\} ). In particular, J^{x,t, k) gives a uniquely defined 
Q{x,t) by 

Q{x,t)= lim -ik[J:3,J+{x,t,k)]. (2.42) 

fc|— )-CXD 

f l2.42p is called the reconstruction formula for Q{x, t). 

Although in general one cannot solve the RH problem explicitely, pure A^-soliton 
solutions can be derived by using the dressing method. We follow the usual approach 
that corresponds to assuming that det a~^{k) (resp. a~{k)) has a finite number N of 
simple zeros kj' G C"*" (resp. k~ G C~), j = 1,...,N. This means that we can use 
the framework of the RH problem with zeroes as described in the previous subsection. 
Considering a trivial regular solution of the RH problem f l2.38p J^{x, t, k) = Ln+i which 
corresponds to Q{x,t) = 0, one then constructs the corresponding RH problem with 
zeroes at k^, j = 1,...,N. The latter enjoys additional properties here due to the 
reduction symmetry (12.50 . Precisely, (I2.17P implies 

kl = {kjy = kj G C+ , (2.43) 

One also gets 

Uj{x, t) = span |e-*<^(-'*''=;)^3 ^f^j^ | ^ ^) = ^^^^ ^) ^ (2.44) 

where Pj is a non-zero vector in C" and Vj~ represents the orthogonal complement of 
Vj. Here Pj is the so-called norming constant associated to kj. Note that fr2.44p implies 
that the projectors involved in the dressing factors are rank-one orthogonal projectors. 
So the degree 1 dressing factor reads 

k — k- 

Di^^{i„„i^^,}{x,t,k) = In+i + {fi,{k) - 1) Ili^^{i,...i^_,}{x,t) , fi^{k) = (2.45) 

and enjoys the property 

^Mn....-.}(^'^'^) = ^l{.i.....-.}(^'^'^*)- (2-46) 

Defining 

0,,{n.....-.}(^,t) = ((^\) , (2.47) 

one gets 

Cij,{n...ij-i}Ci^.,{i^...j^„j}(3^5 ^) 
Ci,-,{n...ij_i}Cij,{n...ij-i}(3^; ^) 

In particular, the reconstruction formula fl2.42p becomes 

N 

Q{x,t) = J2Kkj - A;;)p3,n,-|i,...,,_i|(x,t)] . (2.49) 
i=i 

For later convenience, we introduce the following definition. 



n,,.,...,_,>(a;,t) = . (2.48) 
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Definition 2.9 Let a E Sjy be given and write (cr(l), • • • , cr{N)) = {ii, . . . , ^at). Define 
di^...ii{k), 1 < ^ < A^, recursively by 

, (2.50) 

where, for 1 < j < i, 

di^,{i,...i,_^}{k) = In+ {fi,{k) - 1) 7ri^,{i,...i,^,} , (2.51) 

^ij,{n...ij_i}^*j.{*i---«j-i} 
In the pure iV-soliton system, one has 



TV 



a+{k) = d,„„i^{k) and det a+(A;) = J] /, (fc) , (2.53) 



where a"*" (A;) is defined in (I2.13p . The analytic structure of a"*' (A;) is therefore completed 
determined by the dressing factor di^,,,ij^{k). Finally, we introduce the matrix Aj by 



\im{k-k^){a+{k))-\ deta+(A;,)' - ^^^*«^(^) 



deta+(/cj)' k^kj ^ ' ^ dk 



(2.54) 



Aj contains the information on the residue of {a^{k)) ^ at kj and plays an important 
role for the A^-soliton solution of VNLS on the half- line |l2J . It will be used in section [21 

2.4 Factorization of A^-soliton collision and Yang-Baxter maps 

We are now ready to discuss the factorization property of A^-soliton collisions in the 
VNLS equation. This has already been treated in [H |5] from two different angles. In 
[4j, the factorization is shown by means of an involved direct computation from the 
A^-soliton solution. In [5j, two-soliton collisions are characterized by means of a Yang- 
Baxter map [21]. Although the technical tool is based on asymptotic analyses of the 
A^-soliton solution in both papers and also in ours, we stress that the difference here 
is that our approach is based on the dressing method and its important consequence 
Theorem 12.71 This ensures that the whole construction of the A^-soliton solution is in 
fact consistent from the beginning. In a sense, the factorization of soliton collision is 
a consequence of this consistency. From this point of view, we have an a priori proof 
of factorization whereas the discussions in 11112] can be seen as a posteriori checks that 
the A^-soliton solution obtained with the ISM is consistent. The formulation in terms 
of Yang-Baxter maps turns out to be powerful to discuss factorization and provides a 
natural interpretation of dressing factors and their properties as we will see below. One 
important remark is that strictly speaking, in [5], the map between two polarization 
vectors is established from a two-soliton solution and then factorization is discussed 
based on this. However, to complete the argument, one has to derive such a map 
between two arbitrary polarization vectors within a full A^-soliton solutions. This was 
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in fact the main motivation for the alternative approach of [1]. Here, we perform this 
task directly at the level of the dressing factors using Theorem. 12.71 This results in 
particular in Lemma. 12. 141 below. 

We now make the discussion more precise. First we need to introduce the general 
form of a one-soliton solution characterized by = |(uo + Wq), f o > and /3o, 

^~i{uox+(uQ-vl)t) 

R{x, t) = po vo — — — A — ^^ = Po?o(a;, t) , (2.55) 
cosh(fo(x + 2uot — Axq)) 

where Axq = ^^1^°^ , Po = ]|^- The unit vector po is the polarization of the soliton, 
Wq = —2uq its velocity, vq its amplitude and Axq is the position of the maximum of the 
envelope of the soliton at t = 0. The main feature is that a vector one-soliton is simply 
a polarization vector p times a scalar one-soliton solution q{x,t). 
Now, consider the A^-soliton solution corresponding to 



= o ("j- + ^v,),v,>0 , j = l,...,N, (2.56) 



with the associated norming constants 13 j. Let Wj = —2uj. Then, the following propo- 
sition shows that as t — )■ ±oo, an A^-soliton solution looks like the sum of one-soliton 
solutions up to exponentially vanishing terms. 

Proposition 2.10 Suppose without loss of generality that ui < U2 < ... < mat. De- 
note i?*"(x, t) (resp. i?°"*(x, t) ) the asymptotic solution R{x, t) corresponding tot ^ —oo 
(resp. t ^ oo). Then, 

Rin/out^ t) = y pf"-% ^ + 0(e-^*l*l) . (2.57) 

^ ' ' / J '^1 1 / / J. A in out\\ ^ ' ^ ' 

COSh(fj(X — Wjt — I\Xj )) 

1 tr>in/out\ rtin/out 

jj . ~ . I I A in out ill IP, I 7 in out Pi ..7 

here v = mmt>j, w = mm \ wi — Wj\, Axj = — ana Pj = ^^i„/out^ with 



3 



i-1 N 

= n /^(^.*) dU...Nik,) , /^r = n /^(^;) z?. , (2.58) 

e=i t=j+i 
where ffXk) and di^ i^{k) are defined in Sec. \2.3[ 

Proof: We follow the idea of the scalar case [3] which is based on the evaluation of the 
projectors Ilj^{i j_ij{x, t) as t — t- ±oo. To get the result, it is enough to show that R{x, t) 
approaches the one-soliton solution following the trajectory of a particular soliton I i.e. 
X — Wit = constant, and that it vanishes exponentially for all other directions in the 
(x, t)-plane. But here in the vector case. Theorem. 12.71 is crucial and allows us to write 
the dressing factor in the following form 

Dl,„N = Dl . . . -D^_l,{l...^-2}-D£+l,{l. ..£_!} . . . I JV_;^|-D^ £..Af} ' (2.59) 



11 



where the notation {1 j} means that I is not hsted in {1 ... j}. This means that 

Df^ {1 I N} ^^^^ dressing factor added. Then, recalhng f l2.47[ l2l^8|) . one obtains for 

X — W(t = constant as t — t- — oo 



Di . . . -D^v 1^ I t, k) 



de+i...N{k) 







uT=,m 



(2.60) 



whereas for all other directions, the same calculation yields 0{e '''"1*1). Consequently. 



Ce,{i..I...N}i^y ^) 



-j0(x,i,fc*)S3 / Pi 



UjJimke) 
and the reconstruction formula (I2.49P implies 

C^,{l...i...Af}C|{l..i...Ar}(^'^) 



^1+0(6 



'Vw\t\ 



(2.61) 



Q{x,t) = i{ke - k*^) 



^i{i...i..7v}'^^.{i-i-jv}*^^'^) 

Direct calculation then gives for x — wet = constant as t — — oo 

-i{uiX+{uj-vj)t) 



(2.62) 



R{x,t) = pTve 



cosh(t>^(x — Wit — Ax^")) 



+ 0{e 



(2.63) 



with the various parameters being defined in the proposition. The same technique can 
be applied as t — )■ oo to obtain R°'^^{x,t). m 



Remark 2.11 The order Ui < ■ ■ ■ < means that the relative velocity Wj — Wj+i of 
two consecutive solitons is always positive. Consequently, as t ^ —oo, the solitons are 
distributed along the x-axis in the order 1,2, . . . , N . The picture is reversed as t oo. 
The relative positions of the solitons are therefore completely determined as t ^ ±oo. 

Remark 2.12 Using Theorem. \2. 7\ we could have performed the proof analogously but 
choosing any permutation placing i^ at position N hence giving 

instead of fl2.59p . This corresponds to the possibility that the soliton collisions can occur 
in a different order since we do not know their relative positions at an arbitrary time 
t. However, the final result for would be the same. This is the essence of the 

factorization property. It turns out that this can be made precise by assigning an "in- 
termediate time" polarization vector to each soliton and by considering the effect of a 
two-soliton collision within an A^-soliton solution on the assigned polarization vec- 
tors. The map between the polarization vectors before and after the two-soliton collision 
is a Yang-Baxter map satisfying the set-theoretical Yang-Baxter equation. The mathe- 
matical translation of the factorization property of collisions is therefore an associativity 
property of the operation on polarization vectors given by the Yang-Baxter map. 
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Remark 2.13 The quantity Ax° — AxJ" represents the total position shift incurred by 

soliton j through its collisions with the other solitons. The unit vector p™/""* represents 
the asymptotic polarization vector of soliton j before and after all its collisions with the 
other solitons. From the previous remark, these quantities are independent of the order 
of soliton collisions. We see that by eliminating Pj in (I2.58P that is completely 
determined by through the dressing factors. We note that the obtained relations are 
different from those obtained in the original paper by Manakov in that they do not call 
upon recursively for j < i to obtain The original Manakov 's formula made it 

extremely difficult to see the factorization property and in fact led him to conclude that 
it did not hold. 

To complete the argument and finish the proof of the claims in the previous remarks, 
we define the following "intermediate time" polarization vectors. Let 



7i„{i,+i...i^} = YlM^i) 4+i--iv(^i.)A. 



and 



Pjj,{jj+i...ijv} 



7»j,{ij+i...ijv} 
|7ij,{jj+i..ijv} I 



(2.64) 



(2.65) 



So in particular, p*" = pj^{j_|_i...Ar} and p° = Pj,{i...j-i} and they can be pictorially 
represented as 



Pl,{2...Ar} . . . Pj,{j+l...N} P«,{i+l...Af} ... Pat 



A^-soliton collision 



t —7- — oo 



PAf,{l...Ar-l} • • • Pl,{l...l-1} ■■ Pj,{l...j-1} • • • Pi t — > OO 



X 



We can now formulate the following important lemma. 

Lemma 2.14 Choose kj and ki and assume Uj < Ui. Write for convenience ip = ii . . .iq 
for some q E {1, . . . , N} such that j and I are not in {ii, . . . , iq}. Then 



f*{k*) 

PUiip} = {^n + Uj{K) - l)Pj,{/*p}(Pj,{/ip})^) P/,{v} 



(2.66) 
(2.67) 



where 



*m2 ^ 1 , / (^i ki) |2 



\ki /CjP 



\PMip}\ 
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Proof: From fl2.64p . we have 



pe{i...7V}\{i,ip} pe{i...JV}\{i,/,v} 

7^v}= n/^'(^n 7a,M= WUK) d],^ih)p, 

pG{1...7V}\{/,ip} pG{l...Af}\0,i,ip} 

The relation between d] and c/t . imphes that 

Ip Jlp A 

7;,o V} = /;(^/) (4 + {f*{h) - l)p,-{i^}p]_^,^^) 7/,{v} • 

Introduce 

^ . ^ l7j.{v}l 
l7j,{np}| ' 



we have 



^2 



(2.69) 
(2.70) 



(2.71) 
(2.72) 

(2.73) 



4 + /z(/c;)Pav}p[|i^})7Mnp} 



(2.74) 



It is easy to see that Eij = Eji. Inserting fl2T3|) into flZTTl 1X72]) yields fl2:66l 121671) by 
direct calculation. ■ 



Remark 2.15 The relations defined in Lemma 2.14 have a natural interpretation as an 
"intermediate time" pairwise collision between soliton j and soliton I. Since Wj > Wi 
(uj < Ui), after a certain number of collisions with other solitons (related to the set 
{ip} ), soliton j with polarization Pj,{iip} overtakes soliton I with polarization p«,{ip} and 
acquires polarization Pj,{ip} while soliton I has then polarization p«,{jip}. Pictorially, this 
can be represented as 



Philip} 




P«,{ip} 



Pi,{v} 



X 
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We complete the discussion by rewriting the relations in Lemma. 12.141 in terms of the 
following map acting on CP"^ x CP"'^ onto itself (so that the normalizations in f l2.66p . 
f l2.67p are irrelevant) 

7^(fcl,A:2):(pS^p^^V(pf^P?^), (2-75) 
p(^^) -d + - l) iM!^ p» (2 76) 



p = 4 + 'i^-^ - 1 "V , p . (2.77) 
' V^2-A;i /(pfVpfv 

Then, one finds that the map TZ{ki,k2) is a reversible (parametric) Yang-Baxter map 
i.e. it satisfies 



Uuih, h)ni3{ku k3)n23{k2, ks) = 7^23(A;2, k3)ni3ik^, k3)nu{ki, k2) , (2.78) 

and 

n2iik2,ki)n{ki,k2) = Id. (2.79) 

The rewriting of relations between polarization vectors as a Yang-Baxter map is the basis 
of the argument in [5j. We stress however that the crucial difference of our approach is 
that we obtained the Yang-Baxter map in complete generality for arbitrary polarization 
vectors within a full A^-soliton solution and not just from the two-soliton solution. This 
is similar in spirit to the approach by Tsuchida in [1] but again with the importance 
difference that here, this was made possible by our a priori derivation of Theorem. 12.71 
about dressing factors, instead of an a posteriori derivation from the explicit A^-soliton 
solution. We then recover the following result originally formulated in [3] 

Theorem 2.16 An N-soliton collision in the Manakov model can be factorized into a 
nonlinear superposition of ( ] pairwise collisions in arbitary order. 



3 Factorization with an integrable boundary 

In this section, we use the main results of [2] and combine them with the previous 
construction to discuss factorization of A^-soliton solutions of VNLS on the half-line. 
The idea is that such a solution can be obtained as the restriction to x > of a 2N- 
soliton solution of VNLS on the full line provided that the norming constants and the 
poles kj obey suitable mirror symmetry conditions depending on the boundary conditions 

(fT3|) or dm), (irsj). 

To fix ideas, consider kj = ^{uj + ivj), Vj > 0, j = 1, . . . , 2N and assume that 

Uj > , for j = 1, . . . , iV and Mi < ^2 < ■ ■ ■ < Mat . (3.1) 

^This is introduced properly in section HI 
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This corresponds to the situation where sohtons 1 to are the "real" sohtons (on x > 0) 
and the sohtons + 1 to 2N are the "mirror" sohtons (on x < 0) as t — t- — oo. The 
"real" sohtons have negative velocities so they evolve towards the boundary where they 
meet their "mirror" sohtons which then become the real sohtons. The net result when 
restricted to a; > is that sohtons interact with the boundary and bounce back. This 
looks like 



2A^ , 2A^ - 1 



AT, N 



AT + l 



1,2 



A^, 



t — )■ — oo 



A^ + 1,A^ + 2,...2A^, t-^oo 



where the vertical bar represents the boundary. We can now state the results from [2] 
in a form convenient for our purposes here. 

Proposition 3.1 The N-soliton solution of VNLS on the half-line with boundary con- 
ditions fll.3p or fll.4p . fll.Sp is obtained by the dressing procedure described in Sec. \2.^ 
based on kj and Pj, j = 1,...,2N with kj given by (13. ip for j = 1,...,N and the 
following constraints 



i+N 



where 



M{k*)Aj^ 



M{k) 



-N 



for j 



k — ia 
k + ia 



1,- 



a G 



or 



M{k) 



V 



-1,/ G {l...n}\5^ 



(3.2) 



(3.3) 



(3.4) 



and Aj+N is defined in (I2.54p . 



The crux of the matter is to solve for Pj+n given /3j using (13. 2p . Indeed, the definition 
of Aj+N leads to the following expression 



A 



■j+N 



2N 

n 



ki 



k* 



'An,{1...2N-1} ■ ■ ■ '^j+Ar+l,{l...i+iV}^i+Af,{l...j+Af-l} 



dl 



'j+N-l,{l...j+N-2} 



■ ■ ■ 4 ) (kj+N) 



(3.5) 



where di.^{i^,„i^ 



^}{k) was introduced in Def. 12.91 Hence, the equations in (13. 2p are 
coupled nonlinear equations for the P/s. In Appendix |Al we provide an algorithm 
to solve them. With the knowledge of all the norming constants, we have access to the 
"intermediate-time" polarization vectors which are the basis for the proof of factorization 
on the half-line. They enjoy the following property. 
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Proposition 3.2 Consider 2N polarization vectors as defined in ( [i^. 64\ \2.65\) . Let 

{ii . . An} he a permutation of {1 . . . N} and {ki., ki^+N] Pi-, (3i-+N}je{i...N} be the cor- 
responding "real" and "mirror" poles and norming constants satisfying 



k,^+N = -kl , /3,^/3l^N = M(/c* )A,+iv . 
Then the following relations hold 

Pij+N,{h...iN ii+N...ij^i+N} = m(^ij) Pij,{ij+i...ijv} 5 

with the n X n matrix function m{k) defined by 

' k — ia 



h{k) 



k + ia 



m{k) 



or 



o-i 



n) 



ap = I ,p e S , 

ag = -1 ,q e {I . . .n}\S . 



In particular, 



(3.6) 
(3.7) 

(3.8) 
(3.9) 
(3.10) 



Pj+N,{1...N} = m(A;j)Pj-{i...j...Ar} • 

The proof is long and is given in Appendix [Bl Let us look at these relations from the 
soliton collision viewpoint. As t — )■ — oo, the polarization vectors are ordered as follows 
along the x-axis 



P2Af,{1...2Af-l} 5 P2Ar-l,{1...2Af-2} , ■ ■ ■ 5 PAf+l,{l...Af} 

By applying (13. 7p . this becomes 

m(fcAr) Pat , ... , m(/c2) P2,{3...iv} , m(fci) pi,{2...Af} 



Pl,{2...Af} , P2,{3...Af} , ■ ■ ■ PN , t 



Pi, {2. ..AT} , P2,{3...Af} , ■ ■ ■ PN , t 



-oo 



-oo 



In fact, (13. 7p shows that this picture extends to all "intermediate-time" polarization 
vectors. Therefore, any pairwise collision between soliton ii and ij, described by the 
Yang-Baxter map ^Zi^i.{ki^^ki), is accompanied by a "simultaneous" pairwise collision 
between solitons ij + N and ii + N described by TZi-j^n ii+N{kij+N ■, kii+N), and vice versa. 

Consider now the situation evolving from t — )■ — oo. After a certain number of 
pairwise collisions, soliton j is next to the boundary and the pairwise collision that takes 
place is given by TZj+jy j{kj^iy, kj). The map IZj^N jikj+N.kj) is naturally interpreted 
as the reflection map of soliton j on the boundary. After the reflection, soliton j + N, 
now playing the role of the reflected soliton j, undergoes general pairwise collisions of 
the form IZj+Niikj^N^h) with the remaining "real" solitons and IZq^N j+N{kq+N,kj+N) 
with the "mirror" solitons that travel faster than it. 

We can now define the following reflection map from CP""^ X (C \ zR) to itself that 
describes the change of the polarization vector of one soliton when it interacts with the 
boundary 



i3:(p,A;)^(p,-F) 



(3.11) 
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where 

P = 

Pictorially, this correspond to 

Pj+N,{1...N} 



-j-N} 



(3.12) 



Pj,{j+Ni..:j...N} 




+N j{kj^j\f, kj] 



'j+N,{l...j...N} 




For convenience, we introduce a (generalized) parametric notatioij^ for the reflection 
map: B{k) acting only on vector p. We also define its action on a multiplet of vectors 



Bj{kj) : (pi, . . . , pj, . . . , pjv) t-^ (pi, . . . , p •, . . . , pn) 



(3.13) 



We can now state the main theorem of this paper which guarantees the factorization 
property of the soliton colhsions with a boundary. 



Theorem 3.3 Consider the Yang-Baxter map TZij{ki,kj) defined in \2. 751 -11) and 
Bj{kj) in I[3.11\ - l3.13\} . Then the following set theoretical (parametric) reflection equation 
holds as an identity of maps on 



Bi{k,)'}Z2l{-k;, k,)B2{k2)ni2{k,, k2) = 

7^2l(-fc2^ -kl)B2{k2)n^2{-kl, k2)B^{k{] 



(3.14) 



Moreover, the parametric reflection map B{k) satisfies the following involutive property 

B{-k*)B{k) = Id. (3.15) 

Proof: The involutive property can be found by direct calculation from the definition of 
the refiection map. For the set-theoretical refiection equation, we use the set theoretical 
Yang-Baxter equation together with our mirror image picture. Take kj and kj+2 = 
j = 1,2. Since pairwise collision occurs simultaneously on each side of the boundary, 
there are only 2 possible configurations of collisions and they are identical using the 
Yang-Baxter equation 



"^31(^1, ^3)^32(^3, k2)TZii{ki, A;l)7^42(A;4, ^2)^43(^4, ^3)^12(^:1, ^2) = 

^43(^4, k3)ni2{ki, ^2)7^42(^4, k2)Tl32iki, ^2)7^41(^4, A;i)7^3i(/i;3, ki] 



(3.16) 



^Tliis is defined properly in section |31 



18 



Consider now the following sequence of pairwise collisions 

T^siih, ki)TZ32{h, k2)'R'4i{h, h)'R'42{h, k2)T^Az{h, h)T^i2{h, A;2)(pf\ P?> P?) 
=7^3l(/c3, A;l)7^32(A;3, ^2)^41(^4, ^1)^42(^4, ^2)(pS"\ P2"''? pt'Kpi'^) 
=7^3l(/c3, ki)n32{h, k2)'R'4i{h, A;i)(pS'"\ P2**\ pf'^ pf"'') 

=7^3l(fc3,fcl)(pf^pr^pr^pf^) 

=iP^^\p^^\pl^\p?), (3.17) 

with (p[*\ P2 \ Ps'', P4^) being the initial polarization vectors and (Pi'''', P2^'', Ps^'', pi^^) 
the final polarization vectors. Similarly, consider the following sequence of soliton-soliton 
and soliton-boundary collisions 

i3l(fcl)7^2l(-fc2^ k,)B2{k2)n,2{ku k2){ci?,q^^) 
=Sl(A;l)7^2l(-fc2^fcl)^2(fc2)(qf^qf ) 
=i3l(fcl)7^2l(-A:2^fcO(qf^q^) 
=Si(A:i)(qf\qr^) 

= (qS^\qrV (3.18) 

We claim now that if (qP,q^'^) = (p?,P?), then (qj'^q^^) = (pf ,pf )• From (13T7|) 
and f l3.18p . we have 



,^l*-fc2j (p?)tp« 



' k*-k2 \ q^^^(q?^)t 
,/^t-/^2J (q«)tq? 



pf^^-|/.+ f^4^lC^lpr, (3.21) 



pr^f4.(|i::|ilHiw^pr'. (3.3) 



First, from (13. 19^ I3.20p we have p["^ = q|^**'*, since pj*"* = q^*\ j = 1,2. Then, we 
have p[*"^ = p^"^ and q^"*^ = q^"''. Next, recall that ^4 = —k\ and P4*'' = v(\{k2)P^2^ 
so pf = pf dMD give pS™) = ql™\ Finally, since /tg = -kl and pf ^ = 
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ni{ki)pi^\ relations f l3.23[ I3.24p imply pg^^ = q^^-*. It is easy to check P4'''' = q^^ as 
well. The same argument holds for 

^12(^1, k2)TlA3{h, ^3)^42(^4, ^2)^32(^3, k2)TlAi{h, ki)ni3{ki, A;3)(pf\ p^*\ p^'\ pf^) , 

(3.25) 



and 



(3.26) 



Since iKlT} is equal to IK2^ by IKT^ . then (13778]) is equal to (13726]) . (13J[5]) can be 
showed by the same way by considering the unitary property of Ri2{ki,k2) defined in 
(14. 6 P and the proof is complete. ■ 
The previous proof amounts to the following diagram 



1 2 






— ^ 



1 2 




Note that the Yang-Baxter map ( I2.75P is invariant under the diagonal action of U (n) , 
the set of n X unitary matrices, on CP"^^ x CP"~^ and could therefore be defined on 
CP""^ X CP"~Vf^(^)- However, this is not the case in general for the reflection map 13.111 
except when m is proportional to /„. In this case, the reflection map is proportional 
to Id when acting on the polarization vectors and thus reduces to the identity map 
in CP"'^^. Otherwise, in the case 13.91 with S being a proper subset of {1, . . . , A^}, the 
action of U{n) results in a different reflection map. This is the mathematical translation 
of the physical effects seen in P] on the polarizations when the so-called boundary basis 
does not coincide with the polarization basis. Therefore, we have found two classes of 
reflection maps on CP"^ : the identity map and a family parametrized by U{n) 



Bu{k) : 




t/V[/p, 



(3.27) 



where 



m 



diag{l, 



1,-1 



-1^ 



(3.28) 
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4 Reflection maps 



Let us now present some basic elements of the theory of reflection maps. First we recall 
some definitions of Yang-Baxter maps. Let X be a set and TZ:XxX^XxXa map 
from the Cartesian product of X onto itself. Define TZij : X^ — )■ X^, X^ = X x - ■ - xX, 
as the map acting as TZ on the ith and jth factors of the X-fold cartesian product X^ 
and identically on the others. More precisely, if Tl{x,y) = {f{x,y),g{x,y)), x,y & X, 
then 

^ ^ J 5 '^jj('^l) • • • ? ■^n) ("''I) • • • 5 -^i— 1) /'(•^i? ■^j)) • • • ) di.'^ij ■^i)) "^j+^y • • • ? •^n) ; (^•-^) 
^ ^ J ; '^ij('^l) • • • ; -^n) ("^l; • • • ; -^i— 1; di^'^i^ "^i)' ■ ■ ■ 5 f {.■^iy -^J ) ) "^J + l; • • • ) "^n) • l^'^) 

If 7?.ij satisfy the following Yang-Baxter relation 

7^l27^l37^23 = 7^237^l37^l2 , (4.3) 

then T^jj is called a Yang-Baxter map. In particular, for N = 2, TZi2 = TZ. Let P be the 
permutation map defined by P(x, y) = {y, x) and defined TZ2i = PTZi2P- If, TZ satisfies 

n2in,2 = id, (4.4) 

then TZ is called a reversible Yang-Baxter map. 

It is useful to introduce the so-called parametric Yang-Baxter map TZuiki, k2) which 
is an important special case obtained by considering X x Y instead of X above, where 
Y is the set where ki and k2 live. It satisfies the parametric YB equation 

TZuiki, A;2)7^l3(fcl, ks)n23{k2, k^) = 7^23(A;2, A;3)7^l3(A;l, A;3)7^l2(A;l, /C2) , (4.5) 

and the corresponding reversibility condition reads 

n2i{k2,k^)n{kuk2) = Id. (4.6) 

We have seen that vector soliton collisions in VNLS provide a solution of such a map. 
Before we introduce a general notion of reflection map, we note that (13. lip requires an 
extension of the usual definition of parametric maps. In the general acceptance just 
given, the notation 7?.i2(/i;i, k2) is a short-hand for 

ni2 : {xi, ki] X2, fca) ^ {/{xi, ki; X2, ^2), ^i; g{xi, fci; X2, ^2)5^2) , (4.7) 

meaning that the action of TZ on the set Y is trivial. This is not the case for what we 
defined as a parametric reflection map in (13. lip where we had something of the form: 

B : {x,k) ^ {h{x,k),a{k)) , (4.8) 

where (j{k) = —k* was an involution. Defining S{x,k) = {x,a{k)) from X x Y to 
X xY, we can reconcile the usual notion of parametric map with our context by setting 
B = SB where B is the parametric map B{x, k) = {h{x, k),k). With this definition, the 
parametric refiection equation (I3.14p should read in fact 

5lPl7^2l52P27^l2 = 7^2l52P27^l25l5l , (4.9) 
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as an identity on (X xY) x (X x Y). We also note that the Yang-Baxter map discussed 
in this paper has the property 

^lS27^l2Sl52 = 7^2l . (4.10) 

We now introduce the following general definition. 

Definition 4.1 Given four Yang-Baxter maps TZ^^\ j = 1, 2, 3, 4, a reflection map B is 
a solution of the set-theoretical reflection equation 

B{RfiB2n^^i = n^SB2nf^Bi , (4.ii) 

as an identity on X x X . The reflection map is called involutive if 

BB = Id. (4.12) 

Remark 4.2 // we restrict our attention to involutive reflection maps and reversible 
Yang-Baxter maps, the consistency of the previous definition is ensured by requiring that 
the four Yang-Baxter maps 71^2 are related by TZ^2i — '^21'' — "^12^ — '^12'' ■ This is 
assumed in the rest of the paper. 

Allowing for the above extension of the notion of parametric maps, we can define as 
a special case of the general definition the important class of parametric refiection maps. 

Definition 4.3 Given the parametric Yang-Baxter map TZ^iki, , ki,k2 G Y and an 
involution a : Y ^ Y , a parametric reflection map B{k) is a solution of the parametric 
set-theoretical reflection equation 

i3l(A;l)7^2l(A;l,a(A;2))i32(fc2)7^l2(fcl,fc2) = 7^2l(a(A;l),a(A;2))S2(A;2)7^l2(a(A;l),fc2)Sl(fcl) 

(4.13) 

as an identity on X x X . The reflection map is called involutive if 

B{a{k))B{k) = Id . (4.14) 

Note that in this paper, we explicitely found classes of solution in the case X = CP"^^, 
Y = C* and (T{k) = —k* with 7^12(^1, ^2) being the Yang-Baxter map corresponding to 
VNLS. 

We conclude by defining the notion of transfer maps in analogy with those introduced 
in [21]. Fix N > 2 and define for j = 1, . . . , X the following maps of X^ into itself, 

Tj = TZj+ij . . . TlNjBjTZjN ■ ■ ■ TZjj+iTZjj-i . . . TZjiB'-TZij . . . TZj-ij , (4.15) 
where B^ is a solution of 

Bl7^2li327^l2 = 7^2li327^l2i3l , (4.16) 

and B~ a solution of 

Sl7^l2S27^2l = 7^l2i327^2li3l . (4.17) 

Then one proves by direct (but long) calculation the following result 
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Proposition 4.4 For any reversible Yang-Baxter map TZ, the transfer maps ( I4.15P com- 
mute with each other 

Tj% = %T, , j,£ = l,...,iV. (4.18) 

It is known that the set-theoretical YBE has important connection with the braid 
group (or permutation group in the reversible case) acting on [HE]. Here the formal 
connection of the set-theoretical reflection equation with the finite Artin group (or Weyl 
group in the involutive case) of type BCfq is quite apparent from (I4.16p . It is therefore 
an interesting open question to tackle the construction of general reflection maps in the 
same spirit as the general construction of Yang-Baxter maps in (HUH]. Another interesting 
avenue for further investigation is that of the role of the transfer maps defined above 
from the point of view of Poisson-Lie groups and integrable dynamics along the lines 
described in [22j . 



Appendices 

A Algorithm for the construction of the "mirror" norm- 
ing constants 

We need to solve the coupled nonlinear equations for the defined in (13. 2p . The key 
observation is relying on Theorem. 12.71 Let 

di...2N{k) = did2,{i} . . . d2N,{i...2N-i}{k) . (A.l) 



Recall the definition of Aj+N in (I2.54p . The form of the dressing factor (lA.ip leads to 
the following expression 



Aj+N = II I ^7 I (^"2Af,{1...27V-l} ••• "j+Af+l,{l...i+Af}^J+JV,{l-i+JV-l} 



'^j+Ar_i^{l...j+Ar_2} • • • ^1 ^) (%+Af) ; (^-2) 

where 

dMi-j-i}i^) = In + (^^— ^ - ' (A-3) 

^3,{i-3-i} = 7t / > ^{1....-!} = d\,„,.^{k,)(3j . (A.4) 

Take j = N. Inserting ([Ql [Qj) into implies 

PN{C.2N,{1...2N-l}y = M{k%) TT ( y^^^ ) 7r2N,{1...2N-l} ■ (A. 5) 
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Define a n-complex vector V2Ar by 

^2Ar,{1...2iV-l} 



2N-1 



^2N 



|^2Ar,{1...2Af-l}P 

Combining (1A.5|) and (1A.6|) gives 

^2N,{l...2N-l} 



M{k*^) n 



i=l 



^2N 
W2n\' 



/3 



N ■ 



(A.6) 



(A.7) 



With the knowledge of C,2N,{i...2N-i}, we can thus compute d2N,{i...2N-i}{k). Then take 
j = — 1. From (13. 6p and (lA.2p . we derive 



/5Af-1^27V-l,{1...27V-2} -^(^A^-l) H ( "T^^"^^ 

'^2Af,{1...2iV-l}(^2Af-l) 7r27V_i,{i...2Ar-2} 

Since (i2Ar,{i...2Ar-i}(^) is known, define V2Ar_i by 

^27V-l,{1...2Af-2} 



X 



(A.8) 



V2Ar-l 



|^2Af-l,{1...2Af-2}| 



2N 



,{1...2Af-l} 

ihN-i) n 



i=l,i^2iV-l 



k2N-i - ki 



'2N-1 



k* 



(A.9) 



Combining (lA.Sp and (1A.9P gives 

'C2Af-l,{1...2iV-2} 



V2Af-l 
|v2Ar-lp 



(A. 10) 



We can indeed compute (i2iV-i,{i...2Af-2}(^)- Recursively, taking j = A^ — 2, — 3 up 
to 1, we are able to obtain C,j+N,{i...j+N-i} and (ij+Af,{i...j+Af-i}5 j = 1,...,A^. Since 
{kj, are known, we have the full knowledge of rfj^ji . .j_i}(A;), j = 1, . . . , N as 

well. Therefore, Pj+n can be derived thanks to (lA.4p . 



B Proof of Proposition 13.2 



To prove (13. 7p . we need the mirror symmetry (13. 6p and the permutability property of 
dressing transformations (Theorem. 12. 7p . To avoid tedious notations (the notations are 
already tiresome!), we choose to work with {1 . . . A^} instead of {ii . . An}. We need 
to keep in mind that {l...A^} can be indeed replaced by any permutation of itself 
but indexed by the ordered number from 1 to A^, by taking j — )• ij. Regardless of 
permutations, the relations (13.60 always hold with respect to the indices. Write (ii...2Af 
in the following form 

dl...2N = di . . . dN,{l...N-l}dN+l,{l...N} ■ ■ ■ d2N,{1...2N-l} (B-1) 
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Taking j = in (13. 6p , A2N can be written as 



A2N= n ( fc^^ _ fc* ) ^2jV,{1...2jV-l} C^1.^.2jV-l (^2jv) ■ (B.2) 
i=l ^ * ^ 



Substituting this into fl3.6p gives 

2Af-l - k \ 

/5^^k{i...2iv-i} = ^(^^) n y kZ-k* ) ''2JV,{1...2JV-1} • (B.3) 

Taking the definitions fl2.64[ I2.65P and the mirror symmetry kj^N = —k*, we come to 
the following identification 

72Af,{1...2Ar-l} = ^2Ar,{1...2Af-l} • 

(B.4) 

2N-1 . , ; \ -1 



Define U2Ar by 



\(,2N,{1...2N-1} 



^,{1...2JV-1|'^ 

Inserting (1B.4I - IB.6P into IB. 31 gives 



U27V = 77 = M (fcjv) 77V • (B.6) 



P2N,{1...2N-1} = = m{kN) Pn ■ (B.7) 

U2Af 



Then taking j = — 1, we have 

2N 



A2N-I - Y\ ik'^- k*) '^27V,{1...2Ar-l}^27V-l,{1...27V-2} d^^.2N-2ik2N-l) ■ (B.8) 
i=l,ij^2N-l V 2Af i / 

Substituting this into (13. 6p implies 

d2N,{1...2N-l}ik2N-l)M ^(fcAr-l)/5Ar-lC2Ar-l,{l...2Af-2} = JI ( I~^^nT^ ) ^2A'-l,{1...2Af-2} • 

j=l,j^2iV-l ^ * ^ 

(B.9) 

With (1B.7P and /i;j+Ar = — /i^j, we derive 

d2N,{1...2N-l}{k2N-l)M-\k*^_^) = (^In+ ^^^^ ^ -T*j^ ^ P2N,{1...2N-1}pIn,{1...2N--1}^ M~\k%_^) 

= M-i(A;^_i)4(A;;v-i). (B.IO) 
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According to ( I2.64p . it comes to the following identification 



72Ar-l,{1...2iV-2} = T; -rr~ i2N-l,{l...2N-2} ■ (B.12) 



i^2N-l '^2. 



2Af 



Define 



U2N-1 = ] ^ 7 = M (/i;7v-i)7jv-i,{7V} • (B.13) 

|72Af-l,{1...2Ar-2}| 

Combining (IB. 91 - IB.13P together gives 

P2Ar-l,{1...2Ar-2} = , \ = m(/CAr_i) PAr-l,{Ar} • (B.14) 

Recursively, taking j = — 2, — 3 up to 1, the following relations hold 

V;v,{i...,-i+7v}(w)^"'(^;) = ^"'(^;)<{,+i...7v}(%)' J (B.15) 

Then inserting kj^^ = —k*, 'jj,{j+i...N} and 7j+Ar^{i...j+Ar} into (13. 6p yields 

Pj+N,{j...j-i+N} = m{kj) pj^{j+i..jv} • (B.16) 

By applying the correspondence between {1, . . . , A^} and {ii, . . . , i^r} with j — )■ ij, we 
come to (13. 7p . Equation (13.101 is obtained in the same way by assuming the following 
form of the dressing factor 

di...2N = di. . . dj_i^{i_j_2}dj+i,{i...j-i} ■ ■ ■ '^Ar,{i...3„.Ar_i}'^j,{i...3„.Ar}'^i+Ar,{1...7V} 
• • • dj+N-l,{l...j+N-2}dj+N+l,{l...j+N-l} ■ ■ ■ '^2Af,{l...i+7V...2Af-l}'^j+Af,{l...i+]v...2Af} ' i^-^^) 

This means that dj |^ ^| is the last dressing factor added in the product of the first 
A^ dressing factors, and d^_^^ ^at} is the last dressing factor added in the product 

of the total 2A^ dressing factors. Then applying (13. 7p yields directly (I3.10p . 
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